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ABSTRACT 


The  growth  of  wind  waves  in  eoastal  areas  is  limited  by  the  fetch.  Understanding 
this  sheltering  effect  of  the  coastline  on  the  nearshore  wave  climate  is  of  critical 
importance  for  Navy  operations  (e.g.  amphibious  assault  and  mine  countermeasures)  in 
shallow  water.  Whilst  the  effect  of  fetch  limitation  on  the  development  of  the  wave  field 
is  well  understood,  the  effects  of  bottom  topography,  the  presence  of  swell  and  its 
interaction  with  wind  waves,  the  angle  of  the  wind  relative  to  the  coastline  with  regards 
to  the  change  in  effective  fetch,  and  the  effects  of  atmospheric  stability,  are  not  well 
documented.  This  study  investigates  fetch-limited  wind  wave  growth  by  examining  cases 
where  a  steady  wind  blows  at  various  angles  to  a  straight  coastline,  across  a  continental 
shelf,  in  the  presence  of  swell.  The  observed  wind  wave  growth  for  offshore  winds  is 
consistent  with  previous  observations  in  the  North  Sea.  The  presence  of  energetic  swell 
opposing  the  wind  does  not  have  a  significant  effect  on  the  wind  sea  development. 
Refraction  strongly  affects  the  directional  properties  of  wind  waves  on  the  inner  shelf. 
Observed  wave  growth  rates  agree  well  with  predictions  of  the  WAVEWATCH  III 
operational  wave  prediction  model  forced  with  COAMPS  winds. 
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I.  INTRODUCTION 


The  growth  of  waves  due  to  energy  transfer  from  the  wind  has  been  the  subject  of 
numerous  studies.  Whilst  sophisticated  theories  have  been  developed  for  the  wind- 
induced  growth  and  non-linear  dynamics  of  wind  waves,  very  little  is  known  about  the 
dissipation  of  energy  in  breaking  waves,  and  therefore  the  energy  balance  is  not 
completely  understood  to  this  day.  As  a  result,  forecasting  of  wind  waves  has  relied 
heavily  on  empirical  observations. 

For  a  given  wind  speed,  the  growth  of  wind  waves  is  limited  by  two  variables:  the 
fetch  and  the  duration.  The  fetch  is  the  distance  over  which  the  wind  blows  over  the 
ocean,  and  the  duration  is  the  time  interval  of  steady  wind  conditions.  In  the  open  ocean 
where  the  fetch  is  unlimited,  waves  under  a  steady  wind  continue  to  grow  until  the  energy 
input  to  the  waves  equals  the  energy  dissipated.  At  this  stage,  the  waves  have  reached 
their  maximum  height  for  a  given  wind  speed  and  the  sea  is  considered  to  be  fully 
developed.  On  the  other  hand,  near  a  coast  with  an  offshore  wind,  the  fetch  is  limited  and 
a  fully  developed  sea  does  not  arise.  The  wind  is  unable  to  impart  maximum  energy  to 
the  waves  due  to  this  fetch  limitation  and  the  wave  heights  are  less  than  those  for  a  fully 
developed  sea.  Duration  may  further  limit  wave  growth.  If  a  steady  offshore  wind  blows 
for  a  sufficient  duration,  the  sea  will  reach  a  steady  or  equilibrium  state,  in  which  wave 
heights  increase  with  distance  from  shore,  to  reach  full  development  in  the  open  ocean. 
Although  “ideal”  fetch-limited  conditions  in  which  a  homogenous  wind  blows 
perpendicular  to  a  straight  coastline,  across  a  deep  ocean  in  the  absence  of  swell,  are  rare 
in  nature,  the  topic  has  been  the  subject  of  many  studies  and  has  provided  valuable 
insight  into  the  energy  balance  of  wind  waves.  The  object  of  this  study  is  to  further 
investigate  fetch-limited  wind  wave  growth  on  a  depth-limited  continental  shelf  in  the 
presence  of  swell. 

Phillips  (1958)  recognized  that  waves  cannot  continue  to  grow  indefinitely  under 
the  influence  of  the  wind,  and  hypothesized  the  existence  of  a  saturation  or  equilibrium 
range  of  frequencies  over  which  the  energy  input  and  dissipation  are  balanced  and  the 


1 


waves  no  longer  grow.  Using  a  simple  dimensional  argument,  Phillips  proposed  a  high- 
frequeney  equilibrium  range: 

E{f)  =  ag^f-^  (1) 


where  g  is  gravity  and  a  is  a  numerieal  eonstant.  This  inverse-fifth-power  law  was 
shown  to  be  eonsistent  with  the  observed  roll-off  of  wave  speetra  above  the  speetral  peak. 
However,  rather  than  being  a  universal  eonstant,  the  value  of  the  parameter  a  was  later 
found  to  be  variable,  dependent  on  feteh  (Hasselmann  et  ah,  1973). 


By  fitting  eurves  to  empirieal  data,  Pierson  and  Moskowitz  (1964)  expanded  on 
Phillips’  formula  by  adding  to  the  low  frequeney  end  of  the  speetrum,  obtaining  a 
parametrie  representation  for  a  fully  developed  wind  wave  speetrum: 


E{f)  =  ag\ljtrf-^cxl 
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where  the  peak  frequeney, /„ ,  (eorresponding  to  the  maximum  of  the  speetrum)  is  a 
funetion  of  the  wind  speed.  In  this  parameterisation,  a  steep  forward  faee  for  frequeneies 
less  than  smoothly  matehes  a  high  frequeney  tail  with  slope  of  (-5).  This  widely  used 
parameterisation,  ealled  the  Pierson-Moskowitz  speetrum,  beeame  the  basis  for  further 
developments,  whieh  attempted  to  deseribe  the  growth  of  waves  in  “ideal”  feteh-limited 
eonditions.  Kitaigorodskii  (1962)  proposed  a  similarity  theory  in  whieh  the  evolution  of 
the  non-dimensionalised  wave  energy,  e,  and  peak  frequeney,  v: 
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depend  only  on  the  non-dimensional  feteh: 
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E  is  the  total  sea  surfaee  varianee,  U  is  the  wind  veloeity,  and  x  is  the  feteh. 
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A  comprehensive  investigation  of  the  fetch-limited  growth  of  wind  waves  was 
eondueted  in  the  well-known  Joint  North  Sea  Wave  Project  (JONSWAP)  experiment 
(Hasselmann  et  al.)  in  1973.  The  JONSWAP  experiment  used  data  from  a  eross-shore 
array  of  wave  stations  that  were  loeated  at  varying  fetehes  in  the  North  Sea  to  obtain  a 
parameterisation  for  the  speetral  evolution  of  feteh-limited  wind  waves.  The  so-ealled 
JONSWAP  speetrum  essentially  adds  a  peak  enhaneement  funetion  to  the  Pierson- 
Moskowitz  speetrum: 

/ 

L 

where 

[a,  for/>/^ 

defines  the  width  of  the  spectral  peak  and  y  is  the  ratio  of  the  maximum  spectral  energy 
to  the  maximum  of  the  eorresponding  Pierson-Moskowitz  spectrum.  This  formulation 
contains  five  free  speetral  parameters.  The  peak  frequeney,  and  Phillips  eonstant,  a 
were  sealed  aeeording  to  Kitaigorodskii  (1962),  and  fitted  to  the  JONSWAP 
observations: 

/.  =  3.5x-“^ 

a  =  0.076x'° "" 

For  the  remaining  three  parameters,  eonstant  values  were  proposed: 
y  =  3.3 
a,  =  0.07 
o,  =  0.09 

that  best  fit  the  JONSWAP  observations  (Hasselman  et  ah,  1973).  The  JONSWAP 
parameterisation  was  developed  under  “ideal”  wave  generation  eonditions  where 
stationary  offshore  winds  blow  over  deep  water,  perpendicular  to  a  straight  coastline,  in 
the  absenee  of  swell. 
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The  JONS  WAP  observations  eonfirmed  many  of  the  previously  known  speetral 
charaeteristies  sueh  as  narrow  speetral  peaks,  very  steep  forward  faees  and  an  /“^  high 
frequeney  tail  with  a  feteh-dependent  Phillips  eonstant  (Hasselmann  et  ah,  1973).  The 
spectra  also  exhibited  a  pronounced  overshoot  feature,  which  is  illustrated  in  Figure  1. 
With  increasing  fetch,  a  sharp  spectral  peak  was  observed  to  grow,  while  shifting  towards 
lower  frequency  and  leaving  behind  a  high  frequency  tail.  Energy  at  a  specific  frequency 
increases  to  a  maximum  with  increasing  fetch  and  then  reduces  again  as  the  peak  shifts 
towards  the  lower  end  of  the  spectrum  (Komen  et  ah,  1994).  Non-linear  interactions  were 
proposed  as  the  mechanism  by  which  energy  is  distributed  across  the  frequencies. 

Kahma  (1981)  studied  the  growth  of  the  wave  spectrum  with  fetch  under 
conditions  that  were  similar  to  the  JONS  WAP  experiment  except  that  it  was  conducted  in 
unstable  atmospheric  stratification,  whereas  the  JONSWAP  conditions  were  mainly 
stable.  The  observed  linear  increase  of  the  dimensionless  energy  with  increasing 
dimensionless  fetch  was  similar  to  that  observed  in  JONSWAP,  but  with  energy  levels 
that  were  about  a  factor  of  two  higher  at  all  fetches.  Furthermore,  the  high  frequency  tail 
was  shown  to  better  fit  a  (-4)  power  law  as  opposed  to  the  (-5)  law  in  the  JONSWAP 
spectrum. 

Donelan  et  al  (1985)  conducted  a  comprehensive  study  of  fetch-limited  wave 
growth  on  Lake  Ontario,  Canada.  Observations  were  conducted  over  a  wide  range  of 
conditions  with  near-neutral  atmospheric  stability.  Their  results  supported  Kahma’s 
suggestion  of  a  (-4)  power  law  for  the  high-frequency  tail  of  the  spectrum.  Donelan  et  al. 
proposed  a  slightly  modified  JONSWAP  spectral  model,  incorporating  a  (-4)  tail. 

There  is  still  considerable  debate  as  to  what  characteristic  velocity,  U ,  should  be 
used  in  the  parameterisation  of  wind  wave  spectra.  Kitaigorodskii  (1983)  used  t/„,  the 
wind  speed  at  an  altitude  no  longer  affected  by  the  energy  and  momentum  transfer  from 
the  atmosphere  to  the  waves,  in  the  derivation  of  his  scaling  law.  As  t/„  is  hypothetical, 
it  must  be  replaced  by  some  measured  value  of  wind  speed  in  practice.  The  most  widely 
used  parameter  is  the  friction  velocity,  n*,  which  is  independent  of  height.  However,  in 
the  JONSWAP  experiment,  Hasselmann  et  al.  (1973)  argued  that  the  friction  velocity  is  a 
function  of  fetch  and  thus  should  be  considered  as  an  internal  variable  and  proposed 
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using  the  wind  at  10  m  height  above  the  surfaee,t/io,  as  a  more  appropriate  parameter  for 
scaling.  In  the  present  analysis,  t/jo  was  adopted  as  the  characteristic  wind  velocity. 

Most  studies  of  fetch-limited  wave  growth  concern  a  steady  offshore  wind, 
blowing  perpendicular  to  the  shoreline,  driving  waves  travelling  predominantly  in  the 
wind  direction.  During  the  Lake  Ontario  experiment  (Donelan  et  ah,  1985),  the  effect  of 
the  gradient  of  fetch  about  the  wind  direction  was  investigated  for  different  wind 
conditions.  It  was  found  that  when  the  wind  was  blowing  directly  offshore,  the  direction 
of  the  wind  and  the  waves  agreed.  However,  for  oblique  winds  with  a  large  fetch 
gradient,  disparities  between  the  wind  and  wave  direction  were  as  large  as  50°  and  it  was 
therefore  proposed  that  the  fetch  and  wind  component  in  the  wave  direction  were  more 
suitable  scaling  parameters  instead  of  the  wind  speed  and  fetch  along  the  wind  direction. 

Walsh  et  al.  (1989)  observed  fetch-limited  wave  growth  on  the  U.S.  mid- Atlantic 
coast  using  an  airborne  surface  contour  radar  that  provided  high-resolution  directional 
wave  spectra.  Extensive  measurements  were  taken  to  study  slant  fetch  situations  and  it 
was  found  that  this  situation  produced  an  asymmetrical  directional  wave  spectrum  with  a 
shifted  peak.  The  results  from  this  study  were  in  agreement  with  the  JONSWAP  model 
for  wave  growth  with  respect  to  wave  energy  and  peak  frequency  (Walsh  et  ah,  1989) 
and  also  confirmed  the  Donelan  et  al.  (1985)  relationship  for  the  wave  propagation 
direction  in  the  case  of  slant  fetch  situations. 

To  investigate  the  effect  of  water  of  finite  depth  on  the  growth  of  wind  wave 
spectra.  Young  et  al.  (1996)  conducted  an  experiment  in  a  shallow  lake  of  constant  depth 
2  m.  At  short  fetches,  deep-water  waves  existed  and  the  spectral  growth  curves  were 
similar  to  previous  studies.  However,  at  longer  fetches,  the  energy  and  peak  frequency 
became  depth  limited  and  a  single  growth  curve  no  longer  existed.  Hence  they  proposed 
a  family  of  growth  curves  that  depend  on  the  water  depth  relative  to  the  dominant  surface 
wavelength. 

While  investigations  of  idealized  fetch-limited  conditions  have  yielded  much 
insight  into  the  wave  generation  process,  natural  conditions  are  usually  more  complex 
owing  to  unsteadiness  of  the  wind,  presence  of  swell,  and  the  coastline  geometry  and 
bathymetry,  and  require  a  numerical  wave  prediction  model.  The  WAM  model  was  the 
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first  of  today’s  operational  third  generation  wave  models  (WAMDI  Group,  1988),  whieh 
integrate  the  basie  energy  transport  equation  from  first  prineiples  using  defined  souree 
funetions,  without  any  restrietions  on  the  speetral  shape.  The  JONSWAP  observations  of 
feteh-limited  wave  growth  were  used  to  ealibrate  the  souree  terms  in  the  WAM  model. 
The  WAM  model  ineorporates  finite  depth  effeets  and  thus  can  be  used  as  both  an  open 
ocean  and  coastal  wave  prediction  model. 

Whilst  the  WAM  model  is  a  vast  improvement  on  previous  generation  models,  it 
still  has  some  shortcomings.  Whereas  in  the  JONSWAP  observations,  the  energy  of  the 
wind  waves  grows  linearly  with  increasing  fetch  until  a  fully  developed  state  is  reached 
where  the  energy  remains  constant  (Komen  et  ah,  1994),  in  the  WAM  model,  the 
transition  from  fetch-limited  to  a  fully  developed  state  is  more  gradual  and  the  model 
does  not  actually  achieve  equilibrium.  Furthermore,  the  pronounced  overshoot  of  the 
spectral  peak  development  observed  in  JONSWAP  (Figure  1)  is  not  well  reproduced  by 
the  WAM  model,  owing  to  the  crude  Discrete  Interaction  Approximation  (DIA)  of  the 
non-linear  interaction  source  term.  Other  operational  third  generation  wave  models,  such 
as  WAVEWATCH  (Tolman,  1991)  and  SWAN  (Booij  et  ah,  1999),  are  based  on  a 
similar  formulation  and  treatment  of  the  source  terms  as  used  in  WAM  model. 
WAVEWATCH  incorporated  new  refinements  in  numerical  techniques  as  well  as 
different  source  terms,  while  SWAN  was  designed  for  specific  applications  in  shallow 
water.  However,  the  wave  growth  in  these  models  was  calibrated  using  the  JONSWAP 
observations,  and  thus  wave  generation  predictions  are  similar. 

Numerous  studies  have  provided  further  insight  into  the  complex  process  of  fetch- 

limited  wave  growth  but  the  effects  of  the  bottom  and  its  topography  on  the  wave  growth 

have  been  largely  neglected.  Also,  previous  experiments  have  been  conducted  primarily 

in  protected  marginal  seas  and  lakes,  in  the  absence  of  long  period  swell.  Very  little  is 

known  about  the  effect  of  the  presence  of  swell  on  wind  wave  growth.  In  nature  the 

perfect  situation  rarely  exists  and  the  wave  growth  is  influenced  and  modified  by  a 

number  of  variable  factors  -  fluctuating  winds,  irregular  coastlines,  slant  fetches,  rough 

bottom  topography,  atmospheric  stability  and  the  presence  of  swell.  This  study  further 

investigates  some  of  these  processes  and  is  organised  as  follows.  A  description  of  the 

various  experiments  and  data  sources  is  given  in  Chapter  II.  Chapter  III  describes  the 
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criteria  used  for  ease  study  seleetion  and  the  four  oases  that  were  ohosen  for  the 
investigation  of  fetch-limited  wind  wave  growth.  A  oomparison  between  the 
observations  from  one  particular  case  and  wind  wave  growth  relationships  from  the 
literature  is  presented  in  Chapter  IV.  In  Chapter  V  the  data  are  oompared  with  COAMPS 
wind  predictions  and  WAVEWATCH  III  wave  prediotions,  using  various  souroe  term 
formulations.  Conolusions  are  given  in  Chapter  VI. 
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Figure  1.  Evolution  of  wave  speetra  with  feteh  for  offshore  winds.  The  spectra  are  labelled 
with  the  fetch  in  kilometres,  (from  Flasselmann  et  ah,  1973.) 
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II.  DATA  AND  METHODS 


A.  FIELD  SITE  AND  EXPERIMENTS 

1.  Study  Area 

The  datasets  used  in  this  study  were  colleeted  during  two  separate  field 
experiments,  SandyDuck  and  SHOWEX,  whieh  were  eondueted  offshore  of  the  United 
States  Army  Corps  of  Engineers  Eield  Research  Eacility  (ERE)  near  Duck,  North 
Carolina,  in  1997  and  1999  respectively.  The  instrument  sites  lie  along  an  east/west 
transect  across  the  wide  continental  shelf  of  the  Mid-Atlantic  Bight  region,  midway 
between  Chesapeake  Bay  and  Cape  Hatteras,  on  the  United  States  East  Coast.  Over  this 
region,  the  sandy  shelf  is  narrow  in  the  vicinity  of  Cape  Hatteras,  widens  to 
approximately  100  km  at  the  study  site  and  continues  to  broaden  to  the  north  near 
Chesapeake  Bay.  Water  depths  at  the  instrument  sites  range  from  20  -  200  m.  The 
bathymetry  is  characterised  by  a  gently  sloping  inner  shelf  ('-'1:250),  and  mid  to  outer 
shelf  regions  with  an  even  gentler  slope  (1:2500)  that  are  covered  with  large  ridge-like 
features  with  horizontal  scales  of  0(5  km)  and  vertical  scales  of  0(5  m).  The  offshore 
shelf-break  is  oriented  approximately  north/south  while  the  coastline,  a  long  chain  of 
barrier  islands,  has  a  slightly  oblique  NNW/SSE  orientation.  In  order  to  define  winds 
with  an  offshore  component  for  this  analysis,  the  baseline  coast  was  taken  to  lie  along  the 
North  Carolina  Outer  Banks  on  a  line  340/160,  from  true  north. 

The  prevailing  wind  direction  off  the  North  Carolina  coast  changes  with  the 
seasons.  During  the  summer,  southwesterly  winds  dominate  the  weather  pattern,  whilst 
during  winter,  the  winds  are  predominantly  from  the  north  or  northeast.  The  stretch  of 
coast  from  Chesapeake  Bay  to  Cape  Hatteras  is  exposed  throughout  the  year  to  strong 
winds  and  rough  seas,  caused  by  both  tropical  and  extra-tropical  storms.  The  extra- 
tropical  storm,  referred  to  as  the  “Nor’easter”,  can  occur  at  any  time  of  year  but  the  most 
frequent  and  intense  storms  generally  develop  during  November  to  March.  The  region 
off  Cape  Hatteras,  where  the  Gulf  Stream  is  in  close  proximity  to  the  coast,  is  a  common 
place  for  development  or  explosive  deepening  of  these  storms.  The  tropical  storm,  or  the 
hurricane,  as  it  is  known  in  these  waters,  can  also  develop  in  any  month  but  the  most 

likely  season  is  from  June  to  November.  Whilst  movement  of  these  hurricanes  can  be 
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unpredictable  and  erratic,  they  generally  tend  to  first  move  towards  the  west,  then  turn 
northwestward  and  finally  recurve  towards  the  northeast  (National  Ocean  Service,  2002). 
SHOWEX  was  marked  by  the  occurrence  of  four  separate  hurricanes,  including  the 
devastating  Hurricane  Floyd,  and  a  number  of  Nor’easter  storms.  A  detailed  account  of 
the  conditions  during  SHOWEX  is  given  in  Ardhuin  et  al.  (2003).  Cold  fronts  can  also 
rapidly  traverse  the  area  from  west  to  east  and  are  more  frequent  during  the  winter 
months.  Their  passage  results  in  squally  southerly  winds  that  quickly  shift  west  or 
northwesterly  as  the  front  passes. 

2,  SHOWEX 

The  Office  of  Naval  Research  (ONR)  Shoaling  Waves  Experiment  (SHOWEX) 
was  conducted  from  September  to  December  1999.  An  array  of  six  Datawell  Directional 
Waverider  buoys  (XI -X6)  was  deployed  on  a  cross-shelf  transect  across  the  width  of  the 
North  Carolina  continental  shelf,  on  10-11  September  1999,  as  shown  in  Figure  2.  Buoy 
X5  broke  loose  on  24  October  and  was  subsequently  redeployed  on  28  October  in  the 
same  position.  Buoy  X6  broke  loose  on  21  September,  was  later  redeployed  on  15 
October  in  a  slightly  different  location,  and  was  lost  again  on  5  December.  All  other 
buoys  were  recovered  on  12-13  December  1999  (Table  1). 


Buoy 

Distance  from 
Shore  (km) 

Depth 

(m) 

Latitude 

(North) 

Longitude 

(West) 

Dates  Deployed 

XI 

5.4 

21 

36"  13.73 

75"  42.28 

13  Sep -13  Dec  99 

X2 

12.8 

24 

36"  13.62 

75"  36.77 

11  Sep -12  Dec  99 

X3 

22.2 

26 

36"  12.27 

75"  29.91 

10  Sep -13  Dec  99 

X4 

37.7 

33 

36"  09.50 

75"  19.53 

10  Sep -13  Dec  99 

X5 

58.2 

39 

36"  07.79 

75"  06.05 

13  Sep -24  Oct  99 

28  Oct -13  Dec  99 

X6 

86.62 

197 

36"  04.98 

75"  47.47 

13  Sep -21  Sep  99 

86.57 

193 

36"  04.97 

74"  47.50 

15  Oct  -  5  Dec  99 

Table  1. 


Deployment  data  for  Datawe' 


1  Directional  Waverider  buoys  during  SHOWEX. 


The  Datawell  Directional  Waverider  buoys  are  approximately  one  metre  in 
diameter  and  are  surface  following,  directional  buoys.  Wave  height  spectra  and 
directional  wave  properties  can  be  calculated  from  the  buoys’  measurements  of 
horizontal  and  vertical  orbital  displacements,  at  the  sea  surface.  The  buoys  sample 
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continuously  for  26  minutes  and  then  internally  proeess  the  data  yielding  wave  frequeney 
spectra  in  the  range  0.02-0.56  Hz,  as  well  as  the  mean  propagation  direetion  and 
direetional  spread  as  functions  of  frequency.  The  resulting  speetra,  together  with  the  raw 
time  series,  are  transmitted  every  30  minutes  to  a  reeeiver  on  shore  via  an  HF  radio  link. 
Buoys  X5  and  X6  were  too  far  from  shore  for  radio  transmission,  and  therefore  were 
equipped  with  internal  data  loggers. 

3.  SandyDuck 

Whilst  most  of  the  ease  studies  were  seleeted  from  data  eolleeted  during 
SHOWEX,  one  case  from  an  earlier  experiment  was  also  used  to  eomplete  the  range  of 
wind  direetions  for  the  analysis.  The  SandyDuek  Nearshore  Field  Experiment  was 
eondueted  on  the  inner  eontinental  shelf  and  beaeh  near  Duek,  North  Carolina,  from 
August  to  November  1997.  Extensive  measurements  of  the  evolution  of  frequeney- 
direetional  wave  speetra  were  eolleeted  with  a  500  m  aperture  triangular  array  of  nine  sea 
floor  pressure  sensors  that  was  deployed  5  km  from  the  shore,  in  20  m  depth  (Eigure  3). 
The  spaeing  of  the  pressure  sensors  was  ehosen  to  allow  for  omni-direetional  wave 
measurements  as  well  as  to  aecurately  resolve  sea  waves  in  the  frequency  range  0.04  Hz 
-  0.16  Hz.  In  order  to  provide  measurements  of  the  high  frequency  portion  of  the 
speetrum  (0.16  Hz  -  0.5  Hz)  that  is  attenuated  at  the  sea  floor,  a  Datawell  Directional 
Waverider  Buoy  was  also  deployed  within  the  perimeter  of  the  array. 

B,  ADDITIONAL  WAVE  DATA 

Permanently  deployed  wave-measuring  instruments,  maintained  by  the  U.S. 
Army  Corps  of  Engineers  Eield  Researeh  Eaeility,  provided  additional  direetional  wave 
measurements  in  this  region.  A  15-element  eoherent  array  of  bottom  mounted  pressure 
gauges,  loeated  in  8  m  depth  (8M  in  Eigure  2,  see  also  Eigure  3)  about  900  m  offshore, 
ean  resolve  frequeney-direetional  speetra  up  to  frequeneies  of  approximately  0.3  Hz. 
Additional  nearshore  wave  measurements  were  provided  by  a  Datawell  Direetional 
Waverider  buoy  (WR(FRE)  in  Eigure  2),  positioned  4  km  offshore  in  15  m  depth. 

The  National  Data  Buoy  Center  (NDBC)  provided  additional  wave  data  used  in 
this  study  (Eigures  2  and  4).  Approximately  40  km  to  the  south  of  the  SHOWEX  array, 
an  infrared  laser  wave  gauge  mounted  on  the  Diamond  Shoals  Eighthouse  (DSEN7)  C- 
MAN  station,  provided  wave  frequeney  speetra  in  the  frequeney  range  0.03  -  0.40  Hz. 
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Directional  wave  data  was  obtained  from  the  NDBC  3 -metre  discus  buoy  stations  located 
to  the  north  of  the  instrument  site,  at  Delaware  Bay  (44009)  and  to  the  east  of  Virginia 
Beach  (44014).  These  instruments  output  wave  frequency  spectra  in  the  range  0.03  to 
0.40  Hz,  as  well  as  mean  wave  direction  and  directional  spread  as  functions  of  frequency. 

C.  WIND  DATA 

In  order  to  select  appropriate  fetch-limited  conditions  for  the  case  studies, 
historical  wind  records  from  six  NDBC  stations  on  the  North  Carolina  continental  shelf 
were  examined.  The  locations  of  these  stations  are  indicated  in  Table  2  and  Figure  4.  The 
continuous  winds  data  archive  for  each  of  the  instruments  provided  records  of  ten-minute 
averaged  wind  speed  (in  m/s)  and  direction  (from  which  the  wind  is  coming,  in  degrees 
clockwise  from  true  north).  In  some  cases,  continuous  winds  data  was  not  available  so 
wind  speed  and  direction  were  obtained  from  the  NDBC  standard  meteorological  data 
archive.  Here,  wind  speed  and  direction  were  averaged  over  an  eight-minute  period  for 
buoys,  a  two-minute  period  for  land  stations  and  this  information  was  reported  hourly. 
Since  wind  observations  at  DUCN7  for  November  and  December  1999  were  not 
available  from  the  NDBC  historical  archive,  the  data  reported  in  the  FRF  meteorological 
data  archive  was  used  instead.  These  observations  of  wind  speed  and  direction  were 
vector  averaged  and  reported  every  34  minutes. 


Station 

ID 

Location 

Latitude 

Longitude 

Station  Type 

Anemometer 
Height 
(above  MSL) 

44009 

Delaware  Bay 

38.46  N 

74.70  W 

3  m  discus  buoy 

5.0  m 

CHLV2 

Chesapeake  Light 

36.91  N 

75.71  W 

C-MAN  station 

43.3  m 

44014 

Virginia  Beach 

36.58  N 

74.84  W 

3  m  discus  buoy 

5.0  m 

DUCN7 

Duck  Pier 

36.18  N 

75.75  W 

C-MAN  station 

20.4  m 

DSLN7 

Diamond  Shoals  Lt 

35.15  N 

75.30  W 

C-MAN  station 

46.6  m 

CLKN7 

Cape  Lookout 

34.62  N 

76.52  W 

C-MAN  station 

14.4  m 

Table  2.  NDBC  permanent  wind  and  wave  stations  on  the  North  Carolina/Virginia 

continental  shelf.  DUCN7  is  located  within  500m  of  the  FRF  8m  array  (8M). 


To  account  for  the  varying  heights  of  the  station  anemometers  above  mean  sea 
level,  the  measured  wind  speed  was  converted  to  a  common  reference  level  of  ten  metres 
using  the  log  wind  profile  for  statically  neutral  conditions: 
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(V) 


where  m*  is  the  friction  velocity,  k  =  0.4  is  the  von  Karman  constant  and  z  =  10  m  is  the 
reference  level  (Stull,  1988).  The  aerodynamic  roughness  length,  Zq,  is  obtained  from 
Charnock’s  relationship  for  an  aerodynamically  rough  surface: 

Zo  =  —  (8) 

k 

where  a  =  0.0144  is  the  Charnock  parameter  and  g=9.81  ms‘^  is  the  acceleration  due  to 
gravity  (Charnock,  1955). 

Additionally,  winds  for  a  short  period  were  also  obtained  from  a  MiniMet 
meteorological  buoy,  which  was  deployed  during  SHOWEX.  The  buoy  was  deployed  by 
Dr.  Fred  Dobson  of  Bedford  Institute  of  Oceanography,  Canada,  at  36*’  10.04  N,  75*’ 
20.49.  W  in  about  30  m  depth  (see  Figure  2).  This  instrument  recorded  hourly 
observations  of  standard  meteorological  data  over  a  one-month  period  from  18  November 
to  17  December  1999. 
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Figure  2.  North  Carolina/Virginia  continental  shelf  bathymetry  with  locations  of  SHOWEX 
cross-shelf  transect  and  other  wind  and  wave  measuring  instruments  (courtesy  of  Dr 
Fabrice  Ardhuin). 
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Figure  3.  SandyDuck  instrument  loeations  (from  Borbash,  1998). 
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Figure  4.  North  Carolina/Virginia  coastline  with  locations  of  NDBC  stations  and  the 
SHOWEX  cross-shelf  transect  (courtesy  of  Dr  Fabrice  Ardhuin). 
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III.  CASE  STUDIES 


A.  SELECTION  CRITERIA 

In  order  to  select  appropriate  cases  of  fetch-limited  wind  wave  growth  from  the 
three  month  long  SHOWEX  experiment,  the  wind  measurements  from  the  six  NDBC 
stations  described  in  Chapter  II  were  examined.  The  wind  criteria  used  for  case  selection 
are: 

•  steadiness,  in  both  speed  and  direction  for  at  least  three  hours; 

•  approximately  uniform  over  the  study  region; 

•  active  wave  generation  conditions,  taken  to  be  wind  speeds  >  5  m/s;  and 

•  wind  directions  ranging  from  offshore  to  alongshore. 

Based  on  the  approximate  -20*^  N  orientation  of  the  coastline,  winds  with  directions  from 
160*’  clockwise  to  340*’(from  true  north)  were  considered  to  have  an  offshore  component. 

Winds  for  the  NDBC  station  at  the  Duck  Pier  (DUCN7)  for  the  entire  SHOWEX 
period  are  shown  in  Eigures  5  and  6.  The  horizontal  lines  in  these  figures  bracket  the 
range  of  wind  directions  with  an  offshore  component  that  correspond  to  fetch-limited 
conditions.  The  observations  include  numerous  offshore  wind  events.  However,  offshore 
wind  conditions  that  typically  occur  during  the  passage  of  fronts,  storms  and  hurricanes 
are  usually  highly  variable  with  a  rapidly  veering  wind  direction.  There  were  only  a  few 
time  periods  during  SHOWEX  with  sufficiently  steady  offshore  or  alongshore  winds,  for 
the  purposes  of  this  study.  These  cases  were  chosen  somewhat  subjectively,  using  the 
above  criteria  and  examining  the  winds  from  all  the  NDBC  stations  in  further  detail. 

Once  suitable  wind  criteria  were  met,  the  wave  frequency  and  directional  spectra 
for  each  case  were  examined  in  order  to  ensure  that  the  wave  field  was  also  steady  over 
the  selected  time  period.  Each  case  was  then  narrowed  down  to  a  three  hour  time  period 
for  further  investigation  into  the  wave  growth  with  increasing  fetch. 
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B,  CASES  SELECTED 

1.  SHOWEX  (1999) 

Three  time  periods  during  SHOWEX  were  seleeted  for  analysis  of  feteh-limited 
wind  wave  growth.  The  details  are  given  in  Table  3.  The  direetion  of  the  winds  varies 
from  almost  northerly  winds  for  Case  1,  to  alongshore  north-northwesterly  winds  for 
Case  2  and  approximately  offshore  westerly  winds  for  Case  3. 


Case 

Date 

Time 

Average  Wind 
Direction 

Average  Wind 
Speed  (m/s) 

Remarks 

1 

5  Oct  99 

1800-2100 

352 

9.4 

NTy  Winds 

2 

1  Dec  99 

0600-0900 

343 

14.0 

Alongshore 
NNWTy  Winds 

3 

3  Nov  99 

1330-1630 

272 

9.3 

Offshore  WTy 
Winds 

Table  3.  SHOWEX  eases  ehosen  for  analysis. 


Whilst  these  three  eases  span  a  range  of  different  wind  direetions,  they  all  lie  in 
the  same  quadrant  having  only  northerly  or  westerly  eomponents.  Suitable  eases,  whieh 
fit  the  selection  criteria  but  had  offshore  winds  from  a  southerly  direction,  were  not  found 
in  the  three  months  of  SHOWEX  data. 

2.  SandyDuck  (1997) 

In  order  to  obtain  a  broader  range  of  offshore  wind  cases  for  analysis,  an 
additional  case  (Table  4)  was  selected  from  the  SandyDuck  field  experiment,  which  took 
place  two  years  earlier  in  the  same  geographic  location.  Eigure  7  shows  the  winds  from 
four  NDBC  stations  for  the  selected  twelve-hour  time  period. 


Case 

Date 

Time 

Average  Wind 
Direction 

Average  Wind 
Speed  (m/s) 

Remarks 

4 

17  Aug  97 

2200-0100 

224 

9.2 

Offshore  SWTy 
Winds 

Table  4.  SandyDuck  case  chosen  for  analysis. 


Although  the  wind  direction  in  this  case  is  remarkably  steady  across  the  shelf  for  the 
entire  twelve-hour  period,  the  wind  speed  varies  spatially  and  temporally  over  the  study 
area  with  strong  winds  (up  to  15  m/s)  at  Diamond  Shoals  located  100  km  south  of  the 
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SandyDuck  experiment  site.  This  ease  is  more  eomplex  than  the  others  owing  to  these 
large  spatial  variations  in  wind  speed  and  illustrates,  as  explained  below,  non-local 
effects  on  fetch-limited  wave  generation. 

C.  CASE  ANALYSIS 

1,  Case  1-5  October  1999 

The  wind  observations  from  each  of  the  NDBC  stations  for  Case  1  are  shown  in 
Figure  8.  This  case  exhibits  a  relatively  steady,  northerly  wind,  whose  direction  is 
marginally  outside  the  range  of  fetch-limited  conditions.  This  alongshore  wind  has  a 
slight  onshore  component  and  subsequently  the  fetch  is  close  to  infinite  for  most  of  the 
continental  shelf.  Fetch  limitations  are  expected  only  at  the  nearshore  buoys  and  array 
where  the  wave  propagation  paths  are  affected  by  refraction  over  shallow  bathymetry. 

Energy  spectra,  mean  wave  direction  and  directional  spread,  versus  frequency, 
averaged  over  a  three  hour  period  (Figure  8),  are  shown  in  Figures  9  and  10,  the  former 
showing  the  inner  shelf  instruments  whilst  the  latter  shows  the  outer  buoys  in  the 
SHOWEX  transect.  The  mean  wind  direction,  obtained  by  averaging  the  wind  directions 
from  the  four  NDBC  stations  closest  to  the  transect  (CHEV2,  DSEN7,  44014  and 
DUCN7)  over  the  three  hour  time  period,  is  indicated  by  the  horizontal  dashed  line  in  the 
mean  direction  panel.  Data  from  X6  is  unavailable  for  this  case  study,  as  this  buoy  broke 
loose  on  21  September  and  was  out  of  service  until  15  October. 

The  energy  density  spectra  show  a  broad  maximum  in  energy  with  no  clear 
separation  between  the  wind  wave  and  swell  peaks.  Energy  levels  at  the  peak  frequency 
increase  only  slightly  with  increasing  distance  from  shore  with  the  exception  of 
significantly  higher  levels  observed  at  buoy  X5.  As  X5  is  situated  near  the  shelf  break 
and  waves  from  the  north  would  have  to  travel  over  the  shallow  continental  shelf  before 
reaching  this  buoy,  this  increase  in  energy  level  is  most  likely  caused  by  local  refraction 
effects.  At  frequencies  above  the  peak,  the  spectra  from  all  instruments  collapse  on  top  of 
each  other.  The  relatively  small  variations  between  observed  wave  spectra  suggests  that 
fully  developed  conditions  persisted  across  the  shelf  Only  at  the  shallowest  instruments, 
(WRERE  and  8M)  a  reduction  in  wind  sea  energy  levels  is  noted  that  may  be  the  result  of 
refraction  and  coastal  sheltering  effects.  The  significant  variability  at  all  buoys  at  the 
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lowest  frequencies  is  considered  to  be  noise,  caused  by  inaccuracies  in  the  buoy  response. 
As  the  energy  level  in  this  region  of  the  spectrum  is  so  low,  measurements  of  directional 
properties  are  not  expected  to  be  reliable. 

The  mean  wave  directions  exhibit  a  similar  pattern  at  all  sites  except  for  the 
closest  inshore  instrument  (8M  array)  and  the  furthest  offshore  buoy  (X5).  As  the  wind 
direction  is  nearly  parallel  to  the  depth  contours,  the  predominantly  alongshore- 
propagating  wind  wave  energy  will  be  affected  by  refraction  over  the  shallow  bathymetry 
inshore,  and  propagate  along  curved  trajectories  towards  the  coast.  This  refraction  effect 
can  be  seen  in  the  observed  mean  directions  at  low  frequencies  (Figure  9),  which  are  at  a 
significant  angle  to  the  wind.  At  the  shallowest  8M  array  site,  the  mean  wave  direction 
deviates  by  about  20*^  from  all  other  locations.  Here  refraction  is  felt  across  a  wider  range 
of  frequencies  causing  direction  shifts  as  large  as  90*^.  At  high  frequencies,  the  mean 
wave  direction  tends  towards  the  mean  wind  direction  for  all  of  the  instruments  except 
buoy  X5.  A  possible  explanation  for  this  difference  is  that  the  winds  observed  at  the 
offshore  44014  and  44009  stations  were  consistently  blowing  at  about  a  20*^  angle  relative 
to  the  winds  measured  further  inshore  (Figure  8). 

The  directional  spreading  estimates  indicate  that  the  wave  spectra  are  narrowest  in 

the  vicinity  of  the  spectral  peak  and  broader  at  both  higher  and  lower  frequencies, 

consistent  with  previous  observations  (e.g..  Young  and  Verhagen,  1995).  The  8M  array 

shows  significantly  less  directional  spread  at  low  frequencies  than  the  other  instruments, 

possibly  owing  to  severe  refraction  in  very  shallow  water  close  to  the  coast.  In  addition 

to  refraction,  limited  fetch  for  offshore  propagating  waves  may  also  affect  the  directional 

properties  close  to  shore.  As  the  mean  wind  direction  is  slightly  onshore,  waves 

generated  in  the  wind  direction,  or  on  the  offshore  side  of  this  direction,  will  experience 

an  unlimited  fetch  and  therefore  grow  to  full  development.  On  the  other  hand,  waves 

generated  on  the  inshore  side  of  the  mean  wind  direction  will  be  severely  fetch-limited 

due  to  the  presence  of  the  coast,  and  this  will  result  in  reduced  growth.  This 

asymmetrical  fetch  effect  is  most  noticeable  at  the  8M  array,  which  is  only  about  900  m 

from  the  coast.  Waves  arrive  at  this  instrument  only  from  offshore  as  the  extreme  fetch 

limitation  on  the  coastal  side  will  limit  wave  generation  in  that  region,  resulting  in  a 

reduced  directional  spread.  This  sheltering  effect  may  also  contribute,  in  addition  to 
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refraction,  to  the  reduced  energy  levels  noted  at  the  8M  array  and  to  a  lesser  extent,  at  the 
WRFRF  buoy. 

Figure  1 1  shows  the  two-dimensional  energy  density  spectrum  (averaged  over  the 
same  three-hour  interval)  for  four  locations,  spanning  the  shelf  from  the  shallow  8M 
array  to  buoy  X5  offshore.  The  spectra  are  similar  across  the  transect  exhibiting  a  broad 
uni-modal  peak  that  is  slanted  towards  normal  incidence  (70*^)  at  low  frequencies.  This 
peak  bifurcates  into  a  bimodal  structure  at  high  frequencies,  with  two  peaks  distributed 
approximately  symmetrically  on  either  side  of  the  mean  wind  direction.  This  transition  to 
a  bimodal  spectrum  at  high  frequencies  has  been  noted  in  a  number  of  earlier  studies  (eg. 
Young  &  Verhagen,  1995  and  Wang  &  Hwang,  2001)  and  is  considered  to  be  a  robust 
feature  of  fetch-limited  wind  wave  spectra.  The  cause  of  this  bimodality  is  believed  to  be 
nonlinear  interactions  that  transfer  energy  from  the  dominant  waves  that  travel  in  the 
wind  direction  to  obliquely  propagating  waves  at  higher  and  lower  frequencies  (for 
further  discussion,  see  Young  and  Verhagen,  1995). 

2,  Case  II  -  1  December  1999 

The  NDBC  station  wind  observations  for  Case  II,  as  well  as  the  winds  from  the 
MiniMet  buoy,  which  was  operational  during  this  period,  are  shown  in  Figure  12.  The 
mean  winds  in  this  case  are  strong,  steady  north-northwesterlies.  The  direction  differs  by 
only  about  lO*’  from  that  in  Case  I  but  the  winds  are  slightly  stronger  and  more 
alongshore,  further  restricting  the  fetch  at  the  inshore  stations.  The  winds  are  generally 
uniform  in  direction  across  the  instrument  transect,  but  there  is  some  spatial  variability  in 
the  wind  speed.  Comparisons  of  the  observed  winds  at  DUCN7,  the  MiniMet  buoy 
(which  is  in  the  vicinity  of  X4),  and  44014,  indicate  a  gradual  increase  across  the  shelf 
from  about  10-12  m/s  near  the  coast  to  15-18  m/s  at  the  shelf  break. 

Figures  13  and  14  show  the  energy  spectrum,  mean  wave  direction  and  directional 
spread.  The  spectra  are  generally  similar  at  high  frequencies  whereas  the  spectral  peak 
level  increases  across  the  outer  shelf  between  buoys  X2  and  X6.  This  is  likely  the  result 
of  the  gradual  increase  in  wind  speed  across  the  shelf,  as  opposed  to  a  coastal  effect. 
Evidence  of  coastal  sheltering  effects  is  noted  at  intermediate  frequencies  (0.15  -  0.35 
Hz)  on  the  inner  shelf.  In  this  range,  the  spectra  for  the  offshore  instruments  (Figure  14) 
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collapse,  but  at  the  inner  shelf  instruments  (Figure  13),  there  is  a  systematic  decrease  in 
energy  levels  towards  the  shore.  These  variations  on  the  inner  shelf  are  likely  eaused  by 
the  effeets  of  refraetion  and  asymmetric  fetch,  both  of  which  tend  to  reduce  energy  levels 
and  shift  the  mean  wave  direetion  towards  onshore  propagation,  as  diseussed  for  Case  I. 

The  effects  of  asymmetric  fetch  are  more  dramatie  in  the  mean  wave  direetion 
and  direetional  spreading  estimates  for  the  inner  shelf  instruments  in  Case  II  than  for 
Case  I.  For  Case  II,  the  mean  wind  direetion  is  almost  alongshore,  eausing  the  fetehes  on 
the  shore  side  of  the  wind  to  be  even  more  signifieantly  redueed  than  that  for  Case  I, 
where  the  wind  was  direeted  slightly  onshore.  This  is  particularly  evident  at  intermediate 
frequeneies  at  the  8M  array.  Whereas  in  Case  I,  the  mean  direetion  differs  from  the  other 
instruments  by  about  10-20*’,  here  the  differenees  are  as  large  as  dO*’  at  0.3  Hz.  The 
directional  spread  at  this  array  is  also  signifieantly  less  than  that  for  the  other  instruments. 

The  systematie  large  shift  in  the  mean  wave  direetion  away  from  the  wind  at 
lower  frequeneies,  partieularly  evident  at  the  inner  shelf  instruments  (Figure  13),  is 
caused  primarily  by  refraetion.  Low  frequeney  waves  travelling  at  large  oblique  angles 
to  the  eoast,  are  strongly  refracted  over  the  sloping  bottom  and  turn  towards  the  eoast.  At 
the  lowest  frequencies,  waves  are  travelling  almost  normal  to  the  wind  direetion.  As  the 
wind  foreing  is  expected  to  be  weak  at  this  large  angle,  the  effect  can  only  be  explained 
by  refraction. 

Figure  15  shows  the  two-dimensional  speetra  for  Case  II  at  four  loeations  aeross 
the  SHOWEX  transect.  The  effect  of  refraction  across  the  shelf  is  notieeable  at  all  buoys 
in  the  form  of  a  slanting  ridge  of  energy.  In  partieular  at  the  8M  array  and  at  XI,  the 
maximum  energy  is  eoneentrated  around  the  onshore  direction,  (--070)  nearly 
perpendieular  to  the  wind  direetion.  Similar  to  Case  I,  the  speetra  show  a  transition  from 
a  unimodal  speetrum  at  low  frequeneies  to  a  bimodal  speetrum  at  high  frequeneies. 
However,  in  this  ease,  the  two  peaks  are  not  symmetrieally  distributed  about  the  mean 
wind  direetion,  but  one  of  the  peaks  is  aligned  with  the  wind  direetion.  Wang  and  Hwang 
(2001)  noted  that  when  the  loeal  wind  direction  and  the  mean  wave  direction  were  not 
elosely  aligned,  an  asymmetry  in  the  bimodal  distribution  resulted.  In  this  ease  there  is 
indeed  eonsiderable  direetional  shear  in  the  energetie  part  of  the  spectrum  (Figure  15). 
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The  nonlinear  interaetions  involving  lower  frequeney  swell  eomponents  travelling  at  a 
large  angle  to  the  dominant  wind  waves,  may  eontribute  to  the  asymmetry  of  the 
speetrum  at  high  frequencies. 

The  differences  between  the  spectra  for  Cases  I  and  II  indicate  that  the  evolution 
of  the  wave  spectrum  is  sensitive  to  the  orientation  of  the  wind  with  regards  to  the 
coastline.  There  is  only  a  lO*’  difference  between  the  mean  wind  directions  of  these  two 
cases,  but  in  Case  II  the  inner  shelf  is  sheltered  by  the  coast,  owing  to  the  combined 
effects  of  restricted  fetch  and  refraction. 

3.  Case  III  -  3  November  999 

Case  III,  with  an  offshore  wind,  features  more  severe  fetch-limitations  and 
refraction  effects  than  the  previous  two  cases.  The  wind  observations  at  all  NDBC 
stations  (Figure  16)  show  a  westerly  wind  that  is  relatively  steady  in  direction,  but  with 
some  variability  (±28%  about  the  mean)  in  speed,  both  spatially  and  temporally.  The 
average  wind  speed  over  the  three-hour  case  time  is  about  9.3  m/s.  Winds  to  the  north,  at 
44009,  are  somewhat  stronger  but  this  does  not  significantly  affect  the  generation  of 
waves  by  offshore  winds  at  the  instrument  transect. 

The  energy  density  spectra  (Figures  17,  18)  are  bimodal  with  clearly  separated 
wind  sea  and  swell  peaks.  The  wind  sea  peaks  exhibit  the  classic  fetch-limited  shape  with 
very  steep  forward  faces,  a  high  frequency  tail  evolving  with  increasing  fetch  with 
pronounced  overshoot  of  the  peak  levels.  In  this  case,  the  mean  wind  direction  (-'270*’)  is 
close  to  perpendicular  to  the  coastline  (orientated  340/160),  resulting  in  limited  fetches 
for  all  wave  generation  directions.  However  the  small  20*^  angle  to  shore  normal  is 
significant,  causing  asymmetric  fetches  for  directions  oblique  to  the  mean  wind  direction. 
Waves  arriving  from  south-westerly  angles  will  have  a  shorter  fetch  than  those  arriving 
from  north-westerly  directions.  As  a  result,  maximum  wave  energy  is  generated  in  a 
direction  shifted  slightly  to  the  north  of  the  wind  direction.  These  offshore  propagating 
waves  travel  at  an  angle  to  the  bathymetry  contours  and  thus  will  undergo  refraction. 
Due  to  the  orientation  of  the  mean  wind  direction  and  the  coastline,  waves  will  tend  to  be 
refracted  towards  the  south. 
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The  combined  effects  of  asymmetrical  fetch  and  refraction  are  evident  from  the 
mean  wave  directions,  particularly  at  the  inner  shelf  instrument  locations  (Figure  17).  At 
high  frequencies,  waves  at  all  locations  travel  approximately  in  the  wind  direction  but  at 
intermediate  frequencies,  the  asymmetrical  fetch  effect  becomes  apparent  and  at  lower 
frequencies,  waves  interact  with  the  bottom  and  are  strongly  refracted  towards  the 
shoreline.  For  example,  at  the  WRFRF  location  (15  m  depth),  0.4  Hz  waves  are  in  deep 
water,  so  the  shoaling  bottom  will  have  no  influence  on  ray  paths.  The  waves  however 
are  observed  to  begin  shifting  from  the  mean  wind  direction  towards  the  south  at  this 
frequency  and  hence  the  effect  of  asymmetric  fetch,  which  causes  more  wave  generation 
at  longer  fetches,  must  be  influencing  wave  propagation  at  this  location  and  frequency. 
Below  0.3  Hz,  waves  are  travelling  almost  normal  to  the  mean  wind  direction.  The  wind 
cannot  drive  waves  at  this  large  an  angle.  At  this  frequency,  waves  are  transitioning  from 
deep  to  shallow  water  and  begin  to  be  influenced  by  the  bottom  and  refracted  towards  the 
coast. 

Figure  19  shows  the  two-dimensional  energy  density  spectra  for  four  instruments 
across  the  shelf  Two  distinct  energy  peaks  are  evident  for  buoys  XI,  X3  and  X6,  which 
indicate  a  clear  separation  of  the  low  frequency  swell  and  the  wind  generated  sea.  The 
effects  of  refraction  are  apparent  in  the  dramatic  shift  in  direction  at  the  intermediate 
frequencies  (0.15  -  0.35  Hz),  in  particular  at  buoy  XI.  Here  the  dominant  wave  direction 
varies  by  about  120*^  from  the  mean  wind  direction  at  0.35  Hz  to  an  onshore  direction  at 
0.2  Hz. 

As  this  case  exhibited  nearly  ideal  fetch-limited  wave  growth  conditions  in  the 
presence  of  swell,  it  was  selected  for  further  investigation  and  comparison  against 
established  growth  curves  and  model  predictions  in  Chapters  IV  and  V. 

4,  Case  IV  -  17  August  1997 

Case  IV,  taken  from  the  SandyDuck  experiment,  exhibits  the  most  extreme  effects 
of  refraction  and  asymmetrical  fetches.  In  this  case,  the  mean  wind  speed  varies  spatially 
over  the  shelf  from  relatively  light  wind  ('-'  6-8  m/s)  inshore  at  DUCN7  but  increasing 
offshore  ('-'  8-12  m/s  at  44014)  and  toward  the  south  with  the  strongest  winds  near  Cape 
Hatteras  ('-'  12-15  m/s  at  DSLN7)  (Figure  7).  As  the  mean  wind  direction  is  from  the 
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southwest  (-224*’)  at  a  large  oblique  angle  to  the  eoastline  (oriented  340/160),  the  fetch  is 
asymmetrical  about  the  wind  direction.  Waves  generated  from  directions  on  the  southern 
side  of  the  mean  wind  direction  will  have  considerably  longer  fetches  than  those 
generated  with  a  more  northerly  component.  As  the  winds  and  fetch  are  greatest  to  the 
south  of  the  instrument  locations,  the  dominant  waves  arriving  at  the  8M  and  20M  arrays 
(locations  shown  in  Figure  3)  are  not  locally  forced,  but  generated  by  the  stronger  winds 
a  few  hundred  kilometres  to  the  south.  These  energetic  waves,  initially  travelling 
offshore  in  the  wind  direction,  undergo  severe  refraction  over  the  shallow  continental 
shelf  bathymetry,  resulting  in  a  complete  turning  of  the  waves  toward  onshore 
propagation  angles  at  the  experiment  site. 

The  strong  refraction  is  evident  in  the  frequency  directional  wave  spectra 
observed  at  the  8M  and  20M  arrays  (Figure  20).  Waves  at  the  20M  array  are  travelling 
almost  alongshore,  towards  the  north,  approximately  normal  to  the  mean  wind  direction. 
The  directional  spread  is  very  small  as  expected  for  the  large  refraction  angle.  By  the 
time  the  waves  reach  the  8M  array,  they  have  been  refracted  further  towards  the  shoreline 
and  arrive  at  a  slight  southerly  angle.  The  slanting  ridge  of  energy  apparent  at  the  8M 
array  is  further  evidence  of  the  significant  refraction  that  is  occurring  at  this  location.  At 
higher  frequencies  the  20M  array  shows  evidence  of  a  bimodal  spectrum,  caused  by  the 
two  different  sources  of  wave  energy  arriving  at  the  array.  One  ridge  of  energy  lies  in  the 
wind  direction  and  is  the  result  of  the  locally  generated  wind  sea.  The  more  energetic 
peak  in  the  alongshore  direction  is  caused  by  the  arrival  of  waves  generated  off  Cape 
Hatteras. 
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September  -  October  1999 


Days 


Figure  5.  Wind  speed  and  direction  (from  true  north)  observations  from  NDBC  station 
DUCN7  during  SHOWEX,  September-October  1999.  The  black  horizontal  lines  on  the 
direction  panel  bracket  the  range  for  offshore  winds  and  the  red  vertical  lines  indicate  the 
cases  selected. 
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November  -  December  1999 


Figure  6.  Wind  speed  and  direction  (from  true  north)  observations  from  NDBC  station 
DUCN7  during  SHOWEX,  November-December  1999.  Same  format  at  Figure  5. 
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17  August  1997 
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Figure  7.  NDBC  station  wind  speed  and  direetion  (from  true  north)  observations  over  a 
twelve-hour  period  on  17  August  1997  (SandyDuck,  Case  IV).  The  blaek  horizontal 
lines  on  the  direetion  panel  bracket  the  range  for  offshore  winds  and  the  black  vertical 
lines  indicate  the  selected  time  interval  for  Case  IV. 
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5  October  1999 
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Figure  8.  NDBC  station  wind  speed  and  direetion  (from  true  north)  observations  over  a 
twelve-hour  period  on  5  October  1999  (SHOWEX,  Case  I).  Same  format  as  Figure  7. 
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5  October  1999  18:00-21:00  (EST) 


Frequency  (Hz) 


Figure  9.  Energy  density,  mean  wave  direction  (in  degrees  from  true  north)  and  directional 
spreading  (in  degrees)  versus  frequency  for  Case  I.  Estimates  shown  for  the  inner  shelf 
instruments  are  based  on  a  three-hour  time  period,  indicated  in  Figure  8.  The  blue  dashed 
line  in  the  mean  direction  panel  indicates  the  mean  wind  direction. 
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5  October  1999  18:00-21:00  (EST) 


Figure  10.  Energy  density,  mean  wave  direction  (in  degrees  from  true  north)  and  directional 
spreading  (in  degrees)  versus  frequency  for  Case  1.  Estimates  shown  for  the  outer  shelf 
instruments  are  based  on  a  three-hour  time  period,  indicated  in  Eigure  8.  Same  format  as 
Eigure  9. 
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8IV1  Array 
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Buoy  X5 
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Wave  Direction  (from) 


Figure  11.  Frequency  directional  wave  spectra  observed  at  the  8M  array  and  buoys  XI,  X4 
and  X5  in  Case  I.  The  black  vertical  lines  bracket  the  range  for  offshore  wind  directions 
and  the  white  dashed  line  indicates  the  mean  wind  direction.  As  the  high  frequency 
portion  of  the  spectrum  is  attenuated  at  the  sea  floor,  the  array  spectra  were  truncated  at 
0.32  Hz. 
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1  December  1999 
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Figure  12.  MiniMet  buoy  and  NDBC  station  wind  speed  and  direetion  (from  true  north) 
observations  over  a  twelve-hour  period  on  1  Deeember  1999  (SHOWEX,  Case  11).  Same 
format  as  Figure  7. 
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1  December  1999  06:00-09:00  (EST) 


Figure  13.  Energy  density,  mean  wave  direetion  (in  degrees  from  true  north)  and  directional 
spreading  (in  degrees)  versus  frequency  for  Case  II.  Estimates  shown  for  the  inner  shelf 
instruments  are  based  on  a  three-hour  time  period,  indicated  in  Eigure  12.  Same  format 
as  Eigure  9. 
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Figure  14.  Energy  density,  mean  wave  direetion  (in  degrees  from  true  north)  and  direetional 
spreading  (in  degrees)  versus  frequency  for  Case  II.  Estimates  shown  for  the  outer  shelf 
instruments  are  based  on  a  three-hour  time  period,  indicated  in  Eigure  12.  Same  format 
as  Eigure  9. 
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Figure  15.  Frequency  directional  wave  spectra  observed  at  the  8M  array  and  buoys  XI,  X3 
and  X6  in  Case  II.  Same  format  as  Figure  1 1 . 
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3  November  1999 


Figure  16.  NDBC  station  wind  speed  and  direetion  (from  true  north)  observations  over  a 
twelve-hour  period  on  3  November  1999  (SHOWEX,  Case  111).  Same  format  as  Figure  7. 
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3  November  1999  13:30-16:30  (EST) 


Figure  17.  Energy  density,  mean  wave  direction  (in  degrees  from  true  north)  and  directional 
spreading  (in  degrees)  versus  frequency  for  Case  111.  Estimates  shown  for  the  inner  shelf 
instruments  are  based  on  a  three-hour  time  period,  indicated  in  Eigure  16.  Same  format 
as  Eigure  9. 
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3  November  1999  13:30-16:30  (EST) 


Figure  18.  Energy  density,  mean  wave  direction  (in  degrees  from  true  north)  and  directional 
spreading  (in  degrees)  versus  frequency  for  Case  111.  Estimates  shown  for  the  outer  shelf 
instruments  are  based  on  a  three-hour  time  period,  indicated  in  Eigure  16.  Same  format 
as  Eigure  9. 
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Figure  19.  Frequency  directional  wave  spectra  observed  at  the  8M  array  and  buoys  XI,  X3 
and  X6  in  Case  Ill.  Same  format  as  Figure  1 1 . 
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Figure  20.  Frequency  directional  wave  spectra  observed  at  the  8M  and  20M  arrays  in  Case 
IV.  The  black  vertical  lines  bracket  the  range  for  offshore  wind  directions  and  the  white 
dashed  line  indicates  the  mean  wind  direction.  As  the  high  frequency  portion  of  the 
spectrum  is  attenuated  at  the  sea  floor,  the  array  spectra  were  truncated  at  0.32  Hz  (8M) 
and  0.16  Hz  (20M).  Data  from  the  co-located  directional  Waverider  at  20M  were  used  to 
estimate  the  directional  spectrum  at  frequencies  higher  than  0.16  Hz. 
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IV.  COMPARISON  WITH  GROWTH  CURVES 


As  Case  III  (3  November  1999)  featured  nearly  “ideal”  offshore  wind  conditions 
similar  to  those  detailed  in  the  JONSWAP  study,  the  growth  rates  and  spectral 
characteristics  are  compared  here  with  published  fetch-limited  growth  relationships. 

The  earliest  and  most  well-known  relationships  for  fetch-limited  wind  wave 
growth  are  those  proposed  by  Hasselmann  et  al.  (1973).  These  consist  of  non- 
dimensionalised  growth  curves  for  energy,  (e)  and  peak  frequency,  (v) : 


£  =  1.6x10-^X 

(9) 

o  c  -0.33 

v  =  3.5x 

(10) 

versus  fetch,  (x),  which  best  fitted  the  JONSWAP  observations.  As  the  JONSWAP  data 
was  obtained  during  mostly  stable  atmospheric  stratification,  Kahma  (1981)  suggested 
that  the  growth  relationship  may  vary,  depending  on  whether  the  atmosphere  is  stable  or 
unstable.  In  the  Bothnian  Sea  experiment,  conducted  during  unstable  stratification, 
stronger  growth  was  observed  with  best-fit  curves: 


£  =  3.5  X 10  '' X 

(11) 

o  —0.33 

V  =  31x  . 

(12) 

Kahma  and  Calkoen  (1992)  further  examined  the  growth  rates  in  numerous  datasets  and 
concluded  that  different  growth  curves  existed  for  stable  and  unstable  atmospheric 
stratification.  They  found  that  the  relationship  between  energy  growth  and  fetch  is 
almost  linear  when  the  atmosphere  is  unstable,  but  the  deviation  from  linear  growth  in  a 
stable  atmosphere  is  significant.  They  proposed  the  following  growth  rates  for  stable 
atmospheric  stratification: 


£  =  9.3x10-’/’'' 

(13) 

v  =  1.9x°’' 

(14) 
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and  for  unstable  stratification: 


£  =  5.4x10-’/‘’"  (15) 

v  =  2.3x■“"^  (16) 


In  order  to  compare  the  observations  for  Case  III  with  the  above  growth  curves, 
the  Bulk  Richardson  Number,  in  the  form  used  by  Kahma  and  Calkoen  (1992)  was 
calculated,  as  a  measure  of  atmospheric  stability: 


g(Ja-TJIZ, 

UU/zf 


(17) 


Here  and  are  the  air  and  sea  temperatures  respectively,  z  is  the  height  of  the  wind 
measurement  and  z,  is  the  height  of  the  temperature  measurement.  The  air  and  sea 
temperatures  were  obtained  from  the  standard  meteorological  data  historical  archive  files 
for  NDBC  station  44014  (NDBC,  2003)  and  these  values,  as  well  as  the  calculated  Bulk 
Richardson  number,  are  shown  in  Figure  21.  Over  the  entire  twelve-hour  period,  the  sea 
temperature  was  consistently  O[10*’  C]  warmer  than  the  air  temperature,  resulting  in  a 
statically  unstable  atmosphere. 

To  construct  growth  curves  of  energy  density  and  peak  frequency  versus  fetch, 
the  fetch  for  each  instrument  location  needed  to  be  determined.  The  average  wind  speed 
and  direction  for  this  case  was  calculated  over  the  three-hour  case  time  (1330-1630) 
using  wind  observations  from  the  four  wind  stations  closest  to  the  SHOWEX  transect  at 
DUCN7,  CHLV2,  DSLN7  and  44014.  As  the  coast  is  a  long,  straight  chain  of  barrier 
islands,  in  order  to  determine  the  fetch,  the  baseline  coast  was  approximated  with  a 
straight  line  along  the  seaward  side  of  the  North  Carolina  Outer  Banks,  oriented  340/160 
from  true  north.  The  fetch  was  estimated  to  be  the  distance  along  the  wind  direction 
(-'270*’),  between  each  instrument  location  and  the  intersection  point  on  the  coastline.  The 
fetch  for  each  location  is  shown  in  Table  5. 
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Instrument 

Fetch  (km) 

WRFRF 

3.93 

XI 

6.04 

X2 

14.20 

X3 

23.55 

X4 

37.20 

X5 

56.21 

X6 

82.05 

Table  5.  The  ealculated  feteh  for  Case  III. 

As  the  observations  used  in  this  study  also  include  a  component  of  swell  energy, 
the  spectra  needed  to  be  partitioned  into  separate  sea  and  swell  components.  For  each 
instrument  a  trough  between  the  sea  and  swell  peaks  was  located,  and  the  energy  to  the 
left  of  this  point  was  removed.  In  the  case  of  X2  to  X6,  the  separation  between  the  two 
peaks  was  quite  distinct,  and  the  spectra  were  easily  separated.  For  buoy  XI  and 
WRFRF  however,  the  region  between  the  swell  and  the  wind  wave  peak  was  rather 
broad,  and  it  was  difficult  to  determine  a  distinct  separation  point.  For  these  two  cases  the 
partitioning  was  done  somewhat  arbitrarily.  The  wind  sea  spectra  for  these  two  locations 
may  therefore  contain  a  small  amount  of  residual  swell  energy.  Figure  22  shows  the 
energy  density  spectrum  across  the  entire  range  of  frequencies,  as  well  as  the  remaining 
“wind  sea”  spectrum  that  was  left  after  removal  of  the  swell  energy.  This  part  of  the 
spectrum  was  then  used  in  the  growth  curve  calculations. 

The  fetch,  energy  and  peak  frequency  were  non-dimensionalised  using  Equations 
3-5  (with  a  mean  wind  speed  of  U  =  9.3  m!s).  The  results  are  compared  in  Figure  23 
with  the  JONSWAP,  Kahma  (1981)  and  (unstable)  Kahma  &  Calkoen  (1992) 
relationships.  The  data  appears  to  best  fit  the  relationships  developed  by  Kahma  (1981) 
from  the  Bothnian  Sea  experiment  where  atmospheric  stratification  was  also  unstable. 

The  good  agreement  with  established  growth  curves  is  surprising  considering  the 
presence  of  energetic  swell  with  a  propagation  direction  that  opposes  the  wind  wave 
direction.  These  results  indicate  that  the  presence  of  swell  may  not  have  a  strong  effect  on 
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fetch-limited  wind  wave  growth,  contradicting  results  in  earlier  laboratory  experiments 
(Holthuijsen,  2000;  and  references  therein). 
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3  November  1999 


Time  (EST) 


Figure  21.  Air  and  sea  temperature,  and  estimated  Bulk  Richardson  Number  at  NDBC 
station  44014  for  a  twelve-hour  time  period  on  3  November  1999.  The  black  vertical 
lines  indicate  the  three-hour  time  period  for  Case  111. 
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Energy  Density  Spectra:  3  November  1999 


Frequency  (Hz) 

Figure  22.  Energy  density  spectra  for  Case  III  (Panel  a.)  and  the  truncated  wind  sea  spectra 
(Panel  b.),  after  the  swell  component  has  been  removed. 
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3  November  1999 


Dimensionless  Fetch 


Figure  23. Dimensionless  energy  and  peak  frequeney  versus  dimensionless  feteh  for 
SHOWEX  observations  over  the  three-hour  time  period  of  Case  III.  Previously 
formulated  growth  relationships  from  the  literature  are  also  shown  for  eomparison. 
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V.  COMPARISON  WITH  MODEL  PREDICTIONS 


Third  generation  numerieal  wave  models  are  the  state  of  the  art  tools  for 
predicting  wind  wave  evolution  across  ocean  basins  and  on  smaller  regional  scales.  This 
type  of  model  is  a  vast  improvement  over  first  generation  models,  which  did  not  account 
for  nonlinear  energy  transfers  across  the  spectrum,  and  second  generation  models,  which 
included  a  simple  parameterisation  for  nonlinear  transfers  and  placed  limitations  on  the 
shape  of  the  wind  sea  spectrum  (Tolman,  1991).  Third  generation  models  provide  a  full 
discrete  spectral  resolution  model  for  the  propagation,  generation  and  dissipation  of  two- 
dimensional  wind  wave  spectra  by  integrating  the  basic  energy  transport  equation: 

— +  V,(c.F)  +  V,(c.F)  =  5  (18) 

dt  XX  k  k 

from  first  principles,  using  defined  source  functions  and  without  placing  any  restrictions 

on  the  spectral  shape.  The  energy  density,  F(t,x,k),  is  a  function  of  time,  t,  the 

horizontal  position  vector,  x,  and  the  wavenumber  vector  k.  The  first  term  in  this 
equation  describes  the  local  rate  of  change  and  the  second  and  third  terms  describe 

propagation  effects  in  terms  of  the  divergence  of  energy  fluxes  in  the  physical,  x(x,y), 

and  spectral,  k(f,d),  space  respectively,  where  and  are  the  corresponding  energy 

X  k 

transport  velocities.  The  right  hand  side  of  Equation  18  represents  the  sum  of  the  source 
terms: 

S  =  +  +  (19) 

which  generally  includes  the  effects  of  wind  input,  5,„,  nonlinear  wave-wave 
interactions,  and  dissipation,  Additional  source  terms  such  as  wave-bottom 
interactions  and  bottom  friction,  ,  can  also  be  added  for  shallow  water  applications 
(see  Komen  et  al.  (1994)  for  a  detailed  account  of  the  theory). 

The  WAM  model  (WAMDI  Group,  1988)  was  the  first  operational  third 
generation  model  and  was  calibrated  against  fetch-limited  wave  growth  data  from  the 

JONSWAP  experiment.  Subsequently  developed  third-generation  models  contain  much 
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of  the  same  physics  and  parameterisations  with  refinements  in  numerical  techniques  and 
alternative  wind  input  and  dissipation  source  terms.  Here  we  use  the  WAVEWATCH  III 
model. 

A,  WAVEWATCH  III 

WAVEWATCH  III  is  an  operational  third  generation  ocean  wave  model  that  was 
developed  by  H.  E.  Tolman  and  colleagues  at  the  U.S.  National  Oceanic  and 
Atmospheric  Administration’s  (NOAA)  National  Centers  for  Environmental  Prediction 
(NCEP).  It  is  a  further  development  of  the  earlier  versions  WAVEWATCH  I,  developed 
at  Delft  University  of  Technology,  and  WAVEWATCH  II,  developed  at  NASA  Goddard 
Space  Elight  Centre,  and  includes  numerous  refinements  in  the  structure  of  the  model,  its 
governing  equations,  numerical  methods  and  physical  parameterisations  (Tolman,  2003). 
WAVEWATCH  III  Version  2.22  was  used  for  model-data  comparisons  in  this  study  and 
a  brief  outline  of  the  model’s  source  terms  is  presented  in  the  following  discussion.  A 
complete  reference  for  the  model  physics,  governing  equations,  numerical  schemes  and 
source  terms  can  be  found  in  the  user  manual  and  system  documentation  of 
WAVEWATCH  III  version  2.22  (Tolman,  2002). 

B,  SOURCE  TERMS 

Whilst  the  nonlinear  interaction  and  wind  input  source  terms  in  Equation  19  are 
based  on  sound  physical  theories  and  have  been  studied  extensively,  the  dissipation  term, 
which  accounts  for  the  loss  of  energy  due  to  wave  breaking,  is  poorly  understood  and 
represented  with  simple  empirical  parameterisations.  WAVEWATCH  III  also  includes  an 
additional  fourth  term,  which  accounts  for  bottom  friction  in  shallow  water 
(Tolman,  2002).  Modifications  and  empirical  fine-tuning  of  the  source  terms  has  led  to 
several  alternative  formulations.  In  the  present  study  we  use  the  well-established  source 
terms  of  the  WAM  model  as  well  as  the  more  recent  formulation  proposed  by  Tolman 
and  Chalikov  (1996).  Differences  between  these  formulations  are  summarised  below. 

1,  Nonlinear  Interactions  (5„,) 

To  accurately  calculate  the  nonlinear  interaction  source  term  in  any  wind  wave 
model,  a  six-dimensional  Boltzmann  integral  must  be  computed.  Since  this  evaluation  is 

computationally  prohibitively  expensive  for  use  in  operational  models,  Hasselmann  et  al. 
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(1985)  developed  the  Discrete  Interaction  Approximation  (DIA),  which  parameterises 
this  term  in  terms  of  a  few  interactions.  The  DIA,  first  used  in  the  WAM  model,  was 
adopted  in  the  WAVEWATCH  111  and  many  other  models. 

2.  Wind  Input  ( 5,„)  and  Dissipation  ( 5^,) 

The  energy  balance  between  the  wind  input  and  dissipation  source  term  is 
considered  to  be  important  for  modelling  wave  growth  as  the  two  source  terms  are 
interconnected,  although  they  represent  two  distinct  processes.  Wind  input  and 
dissipation  are  therefore  treated  as  a  collective  source  term  in  the  WAVEWATCH  111 
model  (Tolman,  2002).  Two  different  forms  of  the  combined  source  terms  are  available 
when  running  the  model.  The  WAM-3  combination  of  source  terms,  used  in  cycles  1  to  3 
of  the  WAM  model  contains  a  wind  input  formulation  based  on  experiments  by  Snyder  et 
al.  (1981)  and  a  simple  semi-empirical  dissipation  term  proposed  by  Komen  et  al.  (1984). 
The  JONSWAP  growth  curves  were  used  to  tune  the  dissipation  source  term  so  that  it 
reproduced  previously  observed  fetch-limited  wind  wave  observations  (WAMDl  Group, 
1988). 

Tolman  and  Chalikov  (1996)  presented  an  alternative  parameterisation  based  on 
the  argument  that  the  dissipation  processes  for  frequencies  at  and  below  the  spectral  peak 
are  inherently  different  from  those  occurring  at  higher  frequencies.  Their  dissipation 
source  term  contains  two  distinct  components,  for  these  two  frequency  regimes.  At  the 
low-frequency  end  of  the  spectrum,  dissipation  assumes  a  form  similar  to  the  dissipation 
of  wave  energy  caused  by  turbulence  in  the  oceanic  boundary  layer.  Eor  the  high- 
frequency  region,  a  diagnostic  parameterisation  is  applied,  which  is  modelled  on  the 
power-law  behaviour  of  the  observed  spectral  tail.  The  resulting  source  term  is  a  linear 
combination  of  these  two  components,  and  for  continuity,  a  transition  zone  is  defined 
between  the  high-  and  low-frequency  components  (see  Tolman  and  Chalikov,  1996,  for 
further  details). 

3,  Bottom  Friction  (S^^,) 

The  source  term  describing  dissipation  due  to  bottom  friction  is  parameterised 
according  to  an  empirical  relationship  derived  from  the  JONSWAP  study  (Hasselmann  et 
ah,  1973).  This  formulation  is  used  in  WAVEWATCH  111  and  WAM  (Tolman,  2002). 
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c. 


MODEL  SET-UP 


During  the  SHOWEX  experiment,  COAMPS  wind  fields  with  0.2*’  resolution 
were  obtained  from  the  U.S.  Navy’s  Fleet  Numerical  Meteorology  and  Oceanography 
Centre  (FNMOC)  and  these  data  were  used  as  the  wind  forcing  in  the  WAVFWATCH  III 
model  for  this  study.  A  six-hour  time  step  was  used  for  the  wind  information,  and 
therefore,  the  00  and  12  UTC  data  values  were  taken  from  the  COAMPS  analysis  cycle, 
whilst  the  06  and  18  UTC  values  were  obtained  from  the  forecast  cycle.  As  the  00  and 
06  UTC  data  on  4  November  1999  were  missing,  the  forecast  from  the  12  UTC  cycle  was 
used  instead  to  interpolate  the  winds  at  these  two  times. 

Figure  24  compares  the  NDBC  observations  of  wind  speed  and  direction,  with  the 
COAMPS  predictions  at  the  nearest  grid  points  for  Case  III.  The  COAMPS  data  show 
relatively  good  agreement  with  the  actual  wind  direction  observations,  but  the  wind 
speeds  are  a  little  more  variable.  The  model  field  at  1300  FST  is  shown  in  Figure  25(a). 
Fven  though  there  are  some  discrepancies,  the  COAMPS  winds  do  capture  the  weaker 
inshore  and  stronger  offshore  winds  as  observed  by  the  NDBC  stations.  In  addition  to 
model  runs  forced  with  COAMPS  winds,  model  runs  forced  with  homogeneous  winds, 
using  the  average  observed  wind  speed  and  direction  over  the  three-hour  case  time  (9.3 
m/s  and  270*’),  were  performed  to  examine  the  importance  of  spatial  variability  in  wind 
conditions. 

The  bathymetry  data  used  in  the  model  grid  was  obtained  from  a  variety  of 
sources  including  the  National  Geophysical  Data  Centre  (NGDC)  high-resolution  digital 
bathymetry  data  archive  and  a  number  of  dedicated  surveys  that  were  conducted  during 
the  DUCK94  Nearshore  Processes  Experiment  and  SHOWEX  (see  Herbers  et  al.  (2000) 
and  Ardhuin  et  al.  (2003)  for  details). 

The  WAVFWATCH  III  model  domain  used  in  this  study  covers  the  area  of  the 
North  Carolina/Virginia  continental  shelf  shown  in  Figure  25(b),  with  a  two-minute  grid 
resolution.  The  wave  spectra  were  discretised  into  24  directions  and  30  frequencies.  The 
directions  start  at  true  north  and  encompass  360*’  clockwise  in  increments  of  15*’.  The 
lowest  frequency  is  0.041  Hz  and  they  are  logarithmically  incremented  by  a  factor  of  1.1. 

A  two-minute  time  step  was  used  for  integrating  both  the  propagation  and  source  terms. 
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In  order  to  account  for  the  arrival  of  swell  from  the  open  ocean,  the  model  was  forced  at 
the  offshore  boundary  with  data  from  buoy  X6.  Swell  frequency-directional  spectra  were 
estimated  from  one-hour  buoy  records  using  the  Maximum  Entropy  Method  (Lygre  and 
Krogstad,  1986).  These  data  were  applied  without  any  time  lag  along  the  offshore 
boundary. 

D,  MODEL/DATA  COMPARISON 

Model-data  comparisons  are  presented  here  for  Case  III,  which  exhibited  classic 
fetch-limited  wave  growth  conditions.  A  number  of  model  runs  were  carried  out  in  order 
examine  the  effects  on  wind  wave  growth  of  different  source  terms  and  also  the  presence 
of  swell.  A  total  of  twelve  runs  with  combinations  of  the  following  model  settings  were 
conducted: 

•  Wind  -  COAMPS  versus  homogeneous  winds  or  no  wind; 

•  Atmospheric  Stability  -  stable  versus  unstable; 

•  Source  Term  -  WAM-3  versus  Tolman  and  Chalikov  (T&C); 

•  Bathymetry  -  flat  bottom  versus  actual  bathymetry;  and 

•  Swell  -  pure  wind  sea  versus  inclusion  of  actual  swell. 

Seven  of  these  runs  listed  in  Table  6  are  discussed  in  detail  in  the  following  section. 


Run 

Description 

Wind 

Stability 

Source 

Term 

Swell 

1 

WAM/Homog 

270/9.3  m/s 

Stable 

WAM-3 

No 

2 

WAM/COAMPS 

COAMPS 

Stable 

WAM-3 

No 

3 

T&C/COAMPS 

COAMPS 

Stable 

T&C 

No 

4 

T&C/COAMPS/Unstable 

COAMPS 

Unstable 

T&C 

No 

5 

WAM/COAMPS/SWELL 

COAMPS 

Stable 

WAM-3 

Yes 

6 

T&C/COAMPS/Unstable/SWEEE 

COAMPS 

Unstable 

T&C 

Yes 

7 

WAM/NO  WIND/SWEEE 

None 

Stable 

WAM-3 

Yes 

Table  6.  Description  of  model  runs. 
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All  of  the  above  listed  runs  were  eondueted  using  the  aetual  bathymetry  over  the 
eontinental  shelf.  Some  additional  runs,  not  shown  here,  were  earried  out  using  a  100 
metre  flat  shelf  but  differenees  in  the  predieted  wind  wave  growth  between  the  runs  with 
and  without  bathymetry  were  small,  with  the  exeeption  of  the  nearshore  region  that  is  not 
well  resolved  by  the  model.  In  the  unstable  atmosphere  runs,  a  eorreetion  faetor  (-10°C) 
was  applied  to  refleet  the  aetual  differenee  between  the  measured  air  and  sea 
temperatures. 

In  order  to  separate  the  swell  and  sea  eontributions  of  the  predieted  speetra,  the 
same  eut-off  frequeney  used  for  separating  the  wind  sea  and  swell  eomponents  in  the 
aetual  data  (diseussed  in  Chapter  IV),  was  applied  to  the  model  results  The  following 
model/data  eomparison  diseussion  is  broken  into  two  parts.  The  first  seetion  eompares 
observations  with  model  results  for  only  the  wind  sea  part  of  the  speetrum,  while  the 
seeond  part  diseusses  the  effeet  of  swell  on  predieted  wind  wave  growth. 

1.  Comparison  of  Model  Predictions  with  Observations  and  Empirical 
Growth  Curves 

In  order  to  eompare  the  results  of  the  model  runs  with  previously  published  feteh- 
limited  growth  relationships,  only  the  wind  sea  eomponent  of  the  speetrum  was 
eonsidered.  Figures  26  and  27  show  dimensionless  energy  and  peak  frequeney 
respeetively  as  a  funetion  of  nondimensional  feteh,  for  the  observed  data,  growth  eurves 
in  the  literature  and  model  Runs  1-4.  Results  using  COAMPS  winds  (Run  2)  are 
generally  in  better  agreement  with  observations  than  predietions  based  on  a  homogeneous 
wind  (Run  1).  The  homogeneous  wind  (Run  1)  better  fits  the  observations,  only  at  the 
two  shortest  fetehes  (WRFRF  and  XI),  whilst  Run  2  whieh  uses  the  COAMPS  winds, 
overestimates  the  energy  growth  and  underestimates  the  peak  frequeney.  There  are  two 
possible  eauses  for  this  diserepaney  near  the  eoast.  First,  the  barrier  islands  are  not  well 
resolved  by  the  COAMPS  wind  model  and  the  eoastal  lee  effeet  may  not  be  well 
represented.  Thus  wind  predietions  near  the  shore  may  be  biased  high,  resulting  in  an 
overestimate  of  the  wave  growth  at  short  fetehes.  Spurious  growth  at  the  shortest  fetehes 
may  also  result  from  the  poor  wave  model  resolution  at  the  eoast.  Buoys  WRFRF  and 
XI  are  only  a  few  grid  points  from  the  eoast,  and  henee,  the  model  may  not  resolve  the 
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wave  growth  and  the  refraetion  over  these  short  fetehes.  Higher  resolution  models  are 
needed  to  aeeurately  represent  the  nearshore  wave  evolution. 

Runs  1  and  2  use  the  WAM-3  souree  term,  whieh  is  known  to  flatten  the  growth 
eurve  eompared  to  the  empirieal  JONSWAP  relationship,  eausing  an  overestimation  of 
energy  at  short  fetehes  and  underestimation  of  energy  at  longer  fetehes.  As  the  peak 
frequeney  evolves  as  the  result  of  non-linear  transfers  to  lower  frequeneies,  an  inerease  in 
energy  will  eause  a  deerease  in  peak  frequeney,  and  henee  the  peak  frequeney  will  be 
underestimated.  These  trends  are  eonsistent  with  the  hindeast  results  shown  in  Figures  26 
and  27. 

In  Run  3,  using  the  alternative  T&C  souree  term,  the  energy  is  underestimated 
aeross  the  entire  speetrum.  Compared  with  the  WAM-3  souree  term,  the  errors  at  longer 
fetehes  are  redueed  and  do  not  inerease  with  feteh.  At  longer  fetehes,  the  slope  is 
eomparable  to  that  of  the  Kahma  and  JONSWAP  relations  and  does  not  tend  to  flatten 
out  as  observed  with  the  WAM-3  souree  term.  When  the  T&C  souree  term  is  foreed  with 
COAMPS  winds,  ineluding  a  eorreetion  for  atmospherie  instability  (Run  4),  the  results 
for  the  growth  of  energy  are  elose  to  the  observed  values.  This  run  appears  to  exhibit  the 
most  aecurate  representation  of  the  observations  aeross  the  entire  shelf.  For  the  peak 
frequeney,  the  shortest  fetehes  are  better  modelled  by  the  T&C  souree  term  (Run  4),  but 
at  longer  fetehes  results  for  the  two  souree  terms  (foreed  with  COAMPS  winds)  are 
similar  (Runs  2  and  4). 

2,  The  Effect  of  Swell 

The  winds  offshore  in  Case  III  oppose  the  energetie  swell  arriving  from  the 
Atlantie  Oeean.  Henee  additional  model  runs  were  eondueted  that  inelude  swell,  in  order 
to  investigate  its  effeet  on  the  predieted  growth  of  the  wind  wave  speetrum.  Figures  28 
and  29  eompare  predietions  of  wind  sea  energy  and  peak  frequeney  in  the  presenee  of 
swell  using  the  WAM-3  (Run  5)  and  T&C  (Run  6,  with  unstable  atmosphere)  souree 
terms  respeetively.  Also  ineluded  is  Run  2,  diseussed  earlier,  whieh  uses  the  same 
WAM-3  souree  term  as  Run  5  but  without  swell,  as  a  benehmark  for  wave  growth  in  the 
absenee  of  swell. 
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By  comparing  Runs  2  and  5  (WAM-3),  it  is  noticeable  that  with  the  inelusion  of 
swell  (Run  5)  the  energy  is  signifieantly  overestimated  across  the  entire  wind  sea 
speetrum  and  the  peak  frequeney  is  underestimated.  To  quantify  the  eontribution  of  swell 
in  the  eombined  wind  sea  and  swell  model  results,  an  additional  Run  7  was  condueted 
without  wind.  The  predieted  swell  energy  levels  in  the  wind  sea  band  are  relatively  small 
(eompared  with  Runs  5  and  7  in  Figure  28),  indieating  that  the  inelusion  of  the  swell 
energy  eauses  the  wind  sea  to  grow  unrealistieally.  This  effeet  has  been  observed  before 
(Holthuijsen  et  ah,  2000)  and  ean  be  explained  by  the  way  in  which  the  WAM-3  source 
term  models  the  dissipation  due  to  whitecapping.  Here,  whiteeapping  depends  on  the 
mean  steepness  of  the  waves.  With  the  addition  of  the  low  frequency  swell  into  the 
speetrum,  the  mean  wave  steepness  is  reduced  and  henee  less  dissipation  is  predicted. 
The  resulting  unbalaneed  wind  input  and  dissipation  eauses  inereased  wind  sea  growth 
that  is  not  observed. 

For  Run  6,  the  Tolman  and  Chalikov  source  term  was  used  in  eonjunetion  with 
the  COAMPS  wind  fields  (including  an  instability  correction)  and  the  presenee  of  a 
swell.  This  run  most  aeeurately  represents  the  physieal  situation  and  yields  the  best 
overall  agreement  with  the  observations  with  the  exeeption  of  a  bias  in  the  peak 
frequeney,  whieh  is  underestimated  at  the  shortest  fetehes  and  overestimated  with 
inereasing  feteh.  Whilst  the  WAM-3  souree  term  shows  a  signifieant  differenee  between 
the  model  results  with  and  without  the  presenee  of  swell,  the  T&C  runs  with  (Run  6)  and 
without  (Run  4)  swell,  yield  similar  results  indieating  that  the  addition  of  swell  has  a 
small  effeet  on  the  wind  wave  growth,  which  is  supported  by  the  observations. 

Predietions  of  energy  density,  mean  wave  direetion  and  directional  spread,  as  a 
funetion  of  frequeney,  at  buoy  X2  using  Wam-3  (Run  5)  and  T&C  (Run  6)  souree  terms, 
are  compared  with  observations  in  Figure  30.  Buoy  X2  was  seleeted  for  this  eomparison 
beeause  it  is  elose  to  the  eoast  so  wave  growth  will  be  influenced  by  refraction,  but 
results  are  unlikely  to  be  affected  by  poor  model  grid  resolution.  Whilst  both  model  runs 
aeeurately  refieet  the  position  of  the  swell  peak,  although  slightly  underestimating  the 
energy  level,  the  wind  sea  region  of  the  speetrum  is  not  well  modelled.  The  observed 
speetrum  exhibits  two  distinet  energy  peaks  representing  the  swell  and  wind  sea  with  a 

smaller  peak  in  between,  possibly  eaused  by  refraction  of  wind  waves  elose  to  the  eoast. 
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Run  6  (T&C)  underestimates  the  energy  levels  at  the  intermediate  frequencies  between 
the  dominant  swell  and  sea  peaks,  showing  a  deep  trough  in  the  spectrum  that  is  not 
observed.  In  comparison,  Run  5  (WAM-3),  tends  to  overestimate  the  energy  across  the 
wind  sea  frequencies  but  does  not  exhibit  such  a  large  discrepancy  at  the  energy 
minimum. 

The  mean  wave  direction  observations  at  buoy  X2  show  the  effects  of  refraction 
at  intermediate  frequencies  (0.2  -  0.3  Hz).  Both  model  runs  clearly  show  the  direction  of 
the  low  frequency  onshore  propagating  swell  and  the  high  frequency  offshore  travelling 
wind  waves,  as  observed,  but  the  region  of  the  spectrum  affected  by  refraction  is  not  well 
modelled.  The  T&C  source  term  apparently  under-predicts  the  growth  of  waves  at  large 
angles  to  the  wind  and  thus  under-predicts  the  spectral  levels  of  the  refracted  components 
below  the  peak  frequency.  On  the  other  hand,  the  WAM-3  source  term  over-predicts  this 
growth  and  therefore  over-predicts  the  spectral  energy.  Directional  spreading  is  not  well 
modelled  by  either  source  term  with  WAM-3  over-predicting  and  T&C  under-predicting 
the  spread. 
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Figure  24.  NDBC  station  wind  speed  and  direction  (from  true  north)  observations  over  a 
fourteen-hour  period  on  3  November  1999  (SHOWEX,  Case  III)  are  compared  with 
COAMPS  model  results.  The  black  horizontal  lines  on  the  direction  panel  bracket  the 
range  for  offshore  winds  and  the  black  vertical  lines  indicate  the  selected  time  interval  for 
Case  III. 
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Figure  25.  COAMPS  wind  speed  and  direetion  predietions  at  1300  ESI  (Panel  a)  and 
WAVEWATCH  111  model  domain  with  signifieant  wave  height  predictions  at  1700  EST 
(Panel  b)  on  3  November  1999  (courtesy  of  Dr  Eabrice  Ardhuin). 
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3  November  1999 


Dimensionless  Fetch 


Figure  26.  Growth  of  dimensionless  energy  as  a  funetion  of  dimensionless  feteh  for  model 
Runs  1-4  (SHOWEX,  Case  III).  The  red  line  indieates  the  JONSWAP  growth  relation 
(Equation  9),  the  blue  dashed  line  is  Kahma’s  growth  curve  (Equation  11)  and  the  green 
dashed  line  is  the  unstable  growth  relation  of  Kahma  and  Calkoen  (Equation  15). 
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3  November  1999 


Dimensionless  Fetch 


Figure  27.  Growth  of  dimensionless  peak  frequency  as  a  function  of  dimensionless  fetch  for 
model  Runs  1-4  (SFIOWEX,  Case  III).  The  red  line  indicates  the  JONS  WAP  growth 
relation  (Equation  9),  the  blue  dashed  line  is  Kahma’s  growth  curve  (Equation  11)  and 
the  green  dashed  line  is  the  unstable  growth  relation  of  Kahma  and  Calkoen  (Equation 
15). 
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3  November  1999 


Figure  28.  Growth  of  dimensionless  energy  as  a  function  of  dimensionless  fetch  for  model 
Runs  2  and  5-7  (SHOWEX,  Case  III).  Same  format  at  Figure  26. 
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3  November  1999 


Figure  29.  Growth  of  dimensionless  peak  frequeney  as  a  function  of  dimensionless  fetch  for 
model  Runs  2  and  5-7  (SHOWEX,  Case  III).  Same  format  as  Figure  27. 
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3  November  1999  13:30-16:30  (EST) 


Figure  30.  Energy  density,  mean  wave  direetion  and  directional  spreading  as  a  function  of 
frequency  for  buoy  X2  and  model  Runs  5  and  6. 
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VI.  CONCLUSIONS 


Whilst  much  is  known  about  the  generation  of  wind  waves  in  “ideal”  eonditions, 
natural  eonditions  are  often  more  eomplex  in  eoastal  environments.  External  influenees 
sueh  as  the  orientation  of  the  wind  with  respect  to  the  eoastline,  the  bottom  topography, 
atmospherie  stability  and  the  presenee  of  swell,  introduce  eonsiderable  variability  into  the 
evolution  of  the  wind  wave  field  that  is  not  well  doeumented.  The  aim  of  this  study  was 
therefore  to  gain  a  greater  understanding  of  feteh-limited  wind  wave  generation  on  a 
natural  eontinental  shelf  exposed  to  swell.  This  study  was  eondueted  over  the  North 
CarolinaWirginia  eontinental  shelf,  off  the  east  coast  of  the  U.S.  in  the  vieinity  of  Duek, 
North  Carolina.  Data  were  eolleeted  with  a  eross-shelf  transeet  of  six  Datawell 
Direetional  Waverider  buoys  during  the  Shoaling  Waves  Experiment  (1999)  and  two 
inner-shelf  instrument  arrays  deployed  during  the  SandyDuck  Nearshore  Eield 
Experiment  (1997).  Additional  wind  and  wave  observations  were  obtained  from 
permanent  stations  in  the  area,  maintained  by  the  U.S.  Army  Corps  of  Engineers  Eield 
Researeh  Eacility  and  the  National  Data  Buoy  Center. 

Eour  ease  studies  with  steady  offshore  or  alongshore  winds  (three  from  SHOWEX 
and  one  from  SandyDuek)  were  seleeted  for  analysis  of  wind  wave  growth  over  the 
eontinental  shelf.  The  observations  revealed  that  the  orientation  of  the  mean  wind 
direetion,  with  respeet  to  the  eoastline,  was  of  critical  importance  in  determining  the 
evolution  of  the  wind  wave  field,  especially  close  to  the  eoast.  for  the  two  oases  with 
alongshore  winds,  wave  properties  were  fairly  uniform  aoross  the  shelf,  indioating  that 
the  waves  were  fully  developed,  exoept  at  looations  olose  to  the  eoast.  Here,  reduced 
energy  levels  were  observed,  a  result  of  eoastal  sheltering  and  refraotion  effeots.  With 
alongshore  winds,  waves  generated  in  the  wind  direetion  or  on  the  offshore  side  of  this 
direetion  experienoe  an  unlimited  fetch  and  grow  to  full  development.  However,  waves 
generated  on  the  inshore  side  of  the  mean  wind  direction  are  severely  feteh-limited,  due 
to  the  close  proximity  of  the  eoastline,  and  experienee  redueed  growth.  This 
asymmetrieal  feteh  effeet  was  most  notieeable  on  the  inner  shelf  where  alongshore- 
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propagating  wind  waves  were  also  affected  by  refraction  over  the  shallow  bathymetry 
and  travelled  along  curved  trajectories  towards  the  coast. 

One  case  with  an  almost  shore-normal  offshore  wind  represented  nearly  ideal 
fetch-limited  conditions,  but  included  energetic  swell  directionally  opposing  the  wind 
seas.  This  case  exhibited  the  classic  fetch-limited  spectral  shape  with  pronounced 
overshoot  of  the  wind  wave  peak  levels,  as  observed  during  JONSWAP.  Estimates  of  the 
Bulk  Richardson  Number  indicated  that  the  atmosphere  was  statically  unstable. 
Consistent  with  this  observation,  the  observed  growth  rates  were  closer  to  the  growth 
curves  proposed  by  Kahma  (1981)  for  unstable  atmospheric  conditions,  than  the 
JONSWAP  growth  relations,  which  were  calibrated  for  a  stable  atmosphere.  The  good 
agreement  with  Kahma’s  growth  curves  indicates  that  the  presence  of  swell  does  not 
significantly  affect  the  evolution  of  the  wind  wave  field.  The  effects  of  coastal  sheltering 
on  directional  wave  properties  were  more  severe  in  this  case,  compared  with  the 
alongshore  wind  cases,  which  was  evident  from  the  large  shift  in  mean  wave  directions  at 
low  frequencies.  Since  the  offshore  wind  had  a  small  downcoast  component,  the 
resulting  offshore  travelling  waves  experience  the  combined  effects  of  refraction  and 
fetch  asymmetry.  As  a  result,  the  lower  frequency  waves  travel  preferentially  downcoast 
at  large  angles  to  the  wind. 

One  case  with  spatially  varying  winds  exhibited  the  effects  of  non-local  wind 
wave  generation  on  the  shelf.  Wind  waves  forced  by  strong  fetch-limited  winds  a  few 
hundred  kilometres  to  the  south  of  the  experiment  site,  were  shown  to  undergo  severe 
refraction  over  the  shallow  bathymetry.  This  resulted  in  a  trapping  of  the  waves  along  the 
shelf  and  a  complete  turning  towards  onshore  arrivals  at  the  instrument  arrays. 

One  offshore  wind  case  was  selected  for  further  examination  and  comparison  with 
model  results.  A  series  of  runs,  with  and  without  swell,  were  conducted  using  the 
WAVEWATCH  III  third  generation  model,  forced  with  COAMPS  winds,  using  two 
different  source  term  formulations.  The  widely  used  WAM  cycle  3  parameterisation  and 
the  more  recent  source  term  proposed  by  Tolman  and  Chalikov  (1996)  only  differ  in  their 
wind  input  and  dissipation  formulations.  Comparisons  with  observations  highlight  some 
of  the  known  shortcomings  of  the  WAM-3  source  terms,  that  is,  a  flattening  out  of  the 
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growth  curves  compared  with  the  JONSWAP  relationship  and  spurious  growth  of  wind 
waves  in  the  presence  of  swell.  Overall,  better  agreement  with  the  observations  was 
found  by  using  the  Tolman  and  Chalikov  source  term,  but  this  formulation  tends  to 
under-predict  the  growth  of  waves  at  large  angles  relative  to  the  wind. 
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